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Abstract. 

^^ | The current standard model of cosmology, the ACDM model, is based on the homogeneous 

f ) ■ FLRW solutions of the Einstein equations to which some perturbations are added to account 

for the CMB features and structure formation at large scales. This model fits rather well 

i-C , the observations provided 95% of the energy density budget of the Universe should be of an 

*-™Hi unknown physical nature, i.e. dark matter and dark energy. Now, the aim of a cosmological 

Q ' model is not merely to reproduce the observations, but also to give a physical understanding 

of the Universe we live in. Moreover, even if the assumption of homogeneity seems to be more 
ryj ' or less valid at large scales, it appears to be in contradiction with observations at intermediate 

Cu ' scales (between the scale of non linear structure formation and that where structures virialize). 

This is the reason why, during the last decade, a community of researchers formed whose aim 
was to look for the best way to take into account the influence of the inhomogeneities seen 
^. ] in the Universe and to construct accurate cosmological models which could possibly get rid of 

■^j- ■ the dark components. This task, which is still in its infancy, is currently progressing towards 

promising results. Two types of methods can be found in the literature: spatial averaging of 
OO , scalar quantities and use of exact inhomogeneous solutions of General Relativity. We will give 

Cn ■ here a brief report of the second one. 
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1. Introduction 

The current standard model of cosmology is based on an hypothesis, the Cosmological 
"Principle", stating that the Universe is spatially homogeneous beyond some large scale. 
However, even if the value of this scale is usually claimed to be of the order of 100 Mpc, 
this assumption has never been proved and is still under debate. A mathematical simplification 
describing the Universe by a FLRW model at all scales is therefore commonly used. Up to 1998, 
this model was Einstein-de Sitter (EdS). 

Then, in 1998, it was found that the Supernova la (SN la) luminosity is dimmer than predicted 
by the EdS model [1]. As a consequence this model was ruled out. However, contrary to what 
was widely advocated at the time, and still is, one cannot put forward accelerated expansion at 
this stage. Actually, all current evidence for dark energy is indirect, because the only effect of 
dark energy is on the expansion history of the Universe. 

What is usually done in cosmology is: first assume a cosmological model, then fit it to data 
measuring observables sensitive to the expansion history through the redshift, thus obtain the 
values of the cosmological parameters of the assumed model. 



In an FLRW cosmology, the expansion rate H = a/a is given by 
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It depends on four constants: the fraction of dark energy, £l\, the fraction of curvature, fi^, 
the fraction of matter, including dark matter, SIm-, the fraction of radiation, Q.r, and one 
function, w\, which is the equation of state parameter of the dark energy. Since H 2 (0) = Hq, 
EjOj = 1, removing one constant, fi^ ~ from the WMAP experiment measuring the CMB 
power spectrum and Qr is also small from the CMB (5 x 10~ 5 ). The two other constants Qa and 
rijv/; and the function w\ are rather well constrained by their effect on the expansion history, 
which is manifest through various cosmological observables. The result is the "concordance" 
model with £l\ ~ 0.7, £Im ~ 0.3 (with ~ 0.25 from dark matter) and w = —1. 

Is it a good cosmological model? It seems to work since it reproduces nearly all the 
cosmological data. But 95% is "dark", i.e. not explained. What about a model whose 
components are so fundamentally not understood and could the accelerated expansion be a 
mere mirage? 

It has been shown shortly after the discovery of the Sn la dimming that inhomogeneities in 
our local Universe can in principle mimic an accelerated expansion. A general demonstration 
for models spherically symmetric around the observer has been provided in [2j. Various 
inhomogeneous models with no cosmological constant reproducing the supernova data have 
also been proposed in [21 [H El [5] ■ 

To understand intuitively how such an effect can arise, it is worth turning to the very 
pedagogical model proposed in [5]. This model is composed of two regions. An inner 
homogeneous region with a given density is matched at some redshift to another homogeneous 
region with a higher density. Both regions decelerate, but since the inner void expands faster 
than the outside region because it contains less matter, the observer inside the void, observing 
a source in the outer region, experiences an apparent acceleration of the Universe expansion. 

Three methods have been proposed to take inhomogeneities into account: i) linear 
perturbation theory which is valid when both the curvature and density contrasts remain 
small (this is not the case in the non-linear regime of structure formation and where the 
SNe la are observed), ii) averaging of scalar quantities, promising and currently improving 
iii) exact inhomogeneous solutions which are valid at all scales and are exact perturbations of 
the Friedmann background which they can reproduce as a limit with any precision. This last 
one is the method we will consider in the following. 

2. One patch Lemaitre-Tolman (LT) models with A = 

2.1. The LT solution 

The Lemaitre[6]-Tolman[7] (LT) solution is a spherically symmetric and non-static solution of 
Einstein's equations, with a dust gravitational source. Its metric in comoving coordinates and 
synchronous time gauge, in units where c = 1, is 

ds 2 = dt 2 ^-^rdr 2 _ R 2 (t,r)(d$ 2 + shindy? 2 ), (2) 

1 + 2E(r) 

where the prime denotes derivation with respect to the radial coordinate r. A first integral of 
the field equations gives a dynamical equation for R: 
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R 2 = 2E + — + -R 2 , (3) 



where the dot denotes derivation with respect to the time coordinate t. The mass density in 
energy units, coming from another field equation, is 

2M' 
**G P=WW , (4) 



Integrating Eq.@ gives 
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t-t B (r). (5) 



2E + 2M/R + ±AR 2 

This solution exhibits three integration functions: the energy per unit mass of the particles 
within the comoving shell of radius r, E(r), the gravitational mass of the particles in that 
shell, M(r), and the "bang time", t B (r) meaning that the singularity (bang or crunch) is not 
simultaneous for different shells. 

To reproduce the cosmological data with no dark energy, the cosmological constant, A, is set 
to zero. Then Eq.© can be solved explicitly and the solutions are 

• when E < (elliptic evolution), 

M 
R(t,r) = — —(I- cos V ), (6) 

(-2E) 3 / 2 
1-sini)= l J (t-t B (r)), (7) 

where ry(t, r) is a parameter 

• when E = (parabolic evolution), 

R(t,r) 

• when E > (hyperbolic evolution): 
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-M(t-t B (r)f 



(8) 



M 
R(t,r) = —(cosh V -l), (9) 

(2E) 3 / 2 
sinh77-77= -± (t-t B (r)). (10) 

Note, however, that a given model can exhibit an E(r) function whose sign changes with r 
while going smoothly through zero. Even if, for simplification sake, the parametric solutions are 
commonly used with a fixed sign for E(r), one must be aware of this possible property of the 
LT models. 

Since all the above formulas are covariant under coordinate transformations of the form 
r = fl , ( r )) one °f the three integration functions can be fixed at our convenience by the choice 
of g. Therefore, each LT solution is fully determined by two arbitrary functions of r and a 
coordinate choice. 
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Figure 1. The void for the GBH central observer. Dotted curve = matter density in units 
of the critical density. Dash-dot curve: Qm- Solid curve: Ht, the transverse expansion rate. 
Dashed curve: Hl, the longitudinal expansion rate. 



2.2. Examples of one patch LT models with A = 

This class of solutions has been first used to construct void models in which a local LT under- 
density centered on the observer and a priori parametrized is matched at some redshift to an 
EdS model. We discuss below two examples of such models. 

In [8], the so-called GBH models are specified by their matter content £7^(0 an d their 
expansion rate H{r) and are governed by five free constant parameters. As we said above, they 
are smoothly matched to an EdS universe at large redshift. These models are fitted to a series 
of observations (CMB first peak, Large Scale Structures (LSS), Baryon Acoustic Oscillations 
(BAO), SN la data, Hubble Space Telescope measure of Hq, age of the globular clusters, gas 
fraction in clusters) to constrain their parameters. They suggest the possibility that the Universe 
spherically smoothed around us might be a large (around 2.5 Gpc) void within an EdS universe, 
with no dark energy. This void is depicted in Fig|IJ with the matter content fijv/ and the 
two expansion rates of the LT model, the transverse expansion rate, Ht = R/R and the 
longitudinal expansion rate Hl = R' jR' . However, even if they are able to reproduce most 
of the observations, they exhibit a value of Hq too low as regards the latest measurements which 
give Hq ~ 74 (km/s)/Mpc [9]. This can be seen in Figj2] where Hq ~ 64 (km/s)/Mpc. 

To try to get rid of this drawback a set of five LT void models have been proposed in |10| . 
They have been tested against the full CMB spectrum, BAO, SN and i^o-The inclusion of a 
nonzero overall curvature drastically improves the goodness of fit of these models, bringing 
them very close to that of the ACDM model. And, by varying the density profile, the value 
of Hq has been increased. However, these models have still difficulty exhibiting a value of Hq 
sufficiently high as regards that of [9|. 

In another work [IT], an LT model has been fitted to two sets of observables - the angular 
diameter distance together with the redshift-space mass density and the angular diameter 
distance together with the expansion rate - assumed to have the same form on our past light 
cone as in the ACDM model. The two left arbitrary LT functions have been determined and 
have given the mass distribution in space-time. Without any a priori assumption on its form, 
the current density profile does not exhibit a void but a hump and the observer is located in a 
shallow and wide funnel on top of this hump, as shown in Figj3j Moreover, in this model, the 
value of Hq is that measured, by construction. 




Figure 2. One of the likelihood for the most simple GBH model (ro characterizes the size of 
the void): in yellow with 1-, 2-, 3-<r contours = the likelihood for the combined data set; in blue, 
purple and green respectively with 1- and 2-<r contours = the likelihood for the individual SN 
la, BAO and CMB data sets. The value of Hq is too low, around 64 (km/s)/Mpc. 
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Figure 3. The current density profile of [IT] does not exhibit a void but a hump. The observer 
is located in a shallow and wide funnel on top of this hump. 



However, single patch LT models with a central observer have been criticized on the ground 
of being in contradiction with the Copernican Principle. Now, when using such models in 
cosmology, one should not claim that the observer is at the very centre of any actual spherically 
symmetric universe. Rather, these models must be considered as mathematical simplifications 
or preliminary steps. Here, the inhomogeneities are smoothed over angles around us and only 
their radial component is considered. Note that this is an improvement as regards FLRW models 
where the inhomogeneities are smoothed over the whole space. Since the exact central location of 
the observer is mandatory by construction, these models cannot take anisotropies into account. 
Effects merely designed to test them, such as CMB dipole and low multipoles, kSZ effect, CMB 



spectral distortion, cosmic parallax etc... can in no way be used to test these models [12] , 

3. Swiss-cheese and meat-ball models 

3.1. Example of an LT Swiss-cheese model 

For a more physical representation of our observed Universe, Swiss-cheese models with no dark 
energy have also been proposed. In [13] the model is a lattice of LT bubbles with radius >350 
Mpc in an EdS background. Initially, the void at the center of each hole is dominated by negative 
curvature and a compensating over-density matches smoothly the density and curvature EdS 
values at the border of the hole. Since the voids expand faster than the cheese, the over-dense 
regions contract and become thin shells at the borders of the bubbles while under-dense regions 
turn into emptier voids, eventually occupying most of the volume. An extra very important 
result is that the evolution of the voids, bends the photon paths and affects more the photon 
physics than the geometry of the inhomogeneities. 

3.2. Quasi- spherical Szekeres Swiss-cheese models 

Another exact solution of Einstein's equation has been used to built Swiss-cheese models, the 
quasi-spherical Szekeres (QSS) solution. 

The Szekeres solutions [13] are dust solutions of Einstein's equations with no symmetry at 
all, i.e. no Killing vector. Their metric in comoving coordinates and synchronous time gauge is 

ds 2 = dt 2 _ e 2a dr 2 _ e 2/3 (dx 2 + ^ (n) 

where a and /3 are functions of (t,r,x,y) to be determined from the Einstein equations. The 
/3' = subfamily has not been considered so far in cosmology, since it does not contain the LT 
solution as a spherically symmetric limit. 

The /?' 7^ subfamily can be parametrized so that a set of 5 independent functions of r 
determines each solution. The sign of one of them fixes the geometry of the {t = const.} 3- 
surfaces and the type of evolution (elliptic, parabolic, hyperbolic). The sign of another function 
determines the geometry of the {t = const, r = const} 2-surfaces which can be quasi-spherical, 
quasi-plane or quasi- hyperbolic. Only the quasi-spherical solution has been found useful in 
cosmology, since it is the only one which exhibits the three Friedmann limits (hyperbolic, flat 
and spherical). This solution possesses an origin, not a centre. 

Its metric in comoving coordinates and synchronous time gauge can be written, with the 
parametrization proposed in |15j . 



^ =tf . (»-yV -^ + ^ (12) 

where $ is a function of t and r, k is a function of r and E is a function of r, x and y: 
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with S(r), P(r), Q(r), functions of r. 

Applying the Einstein equations to the metric (|12p and assuming that the energy momentum 
tensor is that of dust, these equations reduce to the following two: 

&= 2 -f-k + \^\ (14) 

__ 2M' -QME'/E 

Kp - &(*>-*E>iEy (15) 



where k = 8irG, A is the cosmological constant and M is an arbitrary function of r. 
As in the LT solution, the bang time function, t B (r), follows from Eq. (jl4p : 

d$ 

t-t B (r). (16) 



-k + 2M/& + ±A$ 2 

Note that Eq. (|16p exhibits the same form as Eq.([5]) of the LT solution. Therefore, for A = 0, 
the parametric solutions are the same, with the Szekeres function k(r) corresponding to the LT 
function —2E(r), i.e. 



when k > (elliptic evolution), 

k 

k' 
~M 



M 
$(t,r) = — (1-cosr?), (17) 

k 3 / 2 
r,-smrj= -—(t-t B (r)), (18) 



when k = (parabolic evolution), 

$(£,r) 
when k < (hyperbolic evolution): 



9 -M(t-t B (r)) 2 



1/3 



(19) 



M 
$(i,r) = — (1-coshry), (20) 

k 

(-k) z / 2 
sinhTy-Ty= ^ (t-t B (r)). (21) 

Remember that, as for the LT case, the sign of k{r) can change in a given model. 

Since all the formulas given so far are covariant under coordinate transformations of the form 
f = g(r), this means that one of the functions k(r), S(r), P(r), Q(r), M(r) or ts{r) can be fixed 
at our convenience by the choice of g. Hence, each QSS solution is fully determined by only five 
functions of r and a coordinate choice. 

To begin with and for simplicity sake, axially symmetric classes of QSS models with A = 
have been used to construct Swiss-cheese models able to reproduce the supernova data [16]. 
Here, the holes are no more spheres, as in LT Swiss-cheeses, but axially symmetric patches and 
the cheese is either an EdS or an open FLRW background. The luminosity distance is computed 
on axially directed null geodesies. Compared to those obtained with corresponding LT Swiss- 
cheeses, the results of these models are quantitatively different but qualitatively comparable. To 
reproduce the observational data, structures of at least 500 Mpc are necessary. 

3.3. Meat ball models 

The meat ball models are the opposite of the Swiss-cheeses. The holes are replaced by balls and 
the cheese by voids which occupy more volume than the over-densities |17j . Most of the photons 
travel through the voids, missing these localized over-densities. The most interesting finding is 
that the lensing probability distribution function is skewed with a mode at demagnified values. 
It is thus possible that these lensing effects, possibly increased by a selection effect, and also 
because they are stronger in an homogeneous background without cosmological constant than 
in a ACDM model, might explain part of the "dark energy". But these lensing effects alone 
cannot explain the full dark energy. However, added to other effects, they can boost LT models 
or make smaller the volume of voids needed in Swiss-cheeses. 



4. Conclusion and prospects 

The first very simple inhomogeneous models used to deal with the dark energy problem and 
proposed very soon after the discovery of the dimming of the luminosity of the supernovae have 
been with time replaced by more and more sophisticated ones. These are able to reproduce to 
a better accuracy not only the supernova observations but nearly all the available cosmological 
data. 

To try to improve the reproduction of our observed local Universe, from FLRW models to 
axially symmetric Szekeres Swiss-cheeses, the number of symmetries has decreased, aiming at 
the removal of any of them. 

It has been claimed that the geometric, dynamical, lensing effects of inhomogeneities act in 
the same direction [17] . Their addition can therefore modify markedly our view of the ACDM 
universe. 

Even if at the end of the day the effects of inhomogeneities should prove insufficient to explain 
out dark energy, the current results obtained with exact models show that their impact on the 
value of the cosmological parameters will have to be considered as non-negligible. 

Now, to be able to take properly into account all the data available to constrain the models, 
and in particular LSS, BAO and CMB, perturbation theories applicable to the LT and Szekeres 
models will have to be developed. 

The final aim is to construct the metric of the cosmos directly from observations. Preliminary 
codes have been developed for LT models [18| . I19| . [20] . but two main problems remain: the choice 
of the scale and to get rid of spherical symmetry. This implies a jump to full numerical Relativity 
which might be the future aim of cosmological investigations. 
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